Roton-Phonon Interactions in Superfluid 4 He 
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High-resolution neutron resonance spin-echo measurements of superfluid 4 He show that the roton 
energy does not have the same temperature dependence as the inverse lifetime. Diagrammatic 
analysis attributes this to the interaction of rotons with thermally excited phonons via both four- 
and three-particle processes, the latter being allowed by the broken gauge symmetry of the Bose 
condensate. The distinct temperature dependence of the roton energy at low temperatures suggests 
that the net roton-phonon interaction is repulsive. 



PACS numbers: 67.25. dt, 78.70.Nx, 03.75.Kk 

Superfluid 4 He is the archetype of a strongly interact- 
ing Bose-condensed system. Its low-temperature prop- 
erties can be understood in terms of weakly interacting 
quasiparticles called phonons and rotons. At zero tem- 
perature, scattering of phonons and rotons is severely 
restricted by kinematical conditions due to particulari- 
ties of their dispersion depicted in Fig. [T] lj , a feature 
which is intimately related to superfluid flow. At fi- 
nite temperatures, rotons acquire a finite lifetime due 
to scattering from thermally populated rotons by four- 
particle processes (4PP), as predicted in a pioneering 
work by Landau and Khalatnikov (LK) [2 j and sub- 
sequently observed in neutron scattering experiments 
[3]. Surprisingly, manifestations of three-particle pro- 
cesses (3PP) have not been reported for roton dynamics. 
Such particle-nonconserving interactions appear due to a 
macroscopic occupation of the lowest energy level and is 
a hallmark of Bose-Einstein condensation [4] . 

In this Letter, we combine high-resolution neutron 
scattering measurements with theoretical analysis to 
show that roton-phonon processes in superfluid 4 He give 
rise to a distinct temperature dependence of the roton 
energy, A(T), at sufficiently low temperatures. We use 
the neutron resonance spin-echo method [51 [B] , which due 
to its extremely stable spin-echo phase allows to measure 
the energy shift of an excitation with high precision si- 
multaneously with the linewidth. The improved resolu- 
tion in our setup makes it possible to study effects that 
are only revealed at temperatures below 1 K, where the 
roton-roton scattering due to 4PPs predicted by LK no 
longer dominates other processes. The main result of our 
work is that the roton linewidth and energy shift in su- 
perfluid 4 He have different temperature dependencies, in 
particular at low temperatures. This is in contrast to the 
LK theory [5] and subsequent theoretical work [7] , where 
the roton lifetime and energy shift are predicted to have 
the same temperature dependence, both simply propor- 
tional to the thermal roton population. A systematic 
evaluation of the leading self-energy diagrams show that 



the experimental results can be understood in terms of 
the interaction of rotons with thermally excited phonons 
via both three- and four-particle processes. This develop- 
ment opens up new perspectives for studies of quasipar- 
ticle interactions in a variety of Bose-condensed systems. 

The measurements were performed on the TRISP ther- 
mal polarized neutron triple-axis spectrometer at the 
FRM-II reactor (Munich), where radio- frequency flip- 
pers are used to encode the energy of each neutron in 
its Larmor precession angle, thereby improving the en- 
ergy resolution by two orders of magnitude. Our setup 
allows to measure energy shifts with an unprecedented 
precision of 2 mK (~ 0.2 fieV) and lincwidths down to 
1.5 mK (~ 0.17 /xeV). The polarization and energy of 
the incoming neutrons were determined by a supermir- 
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FIG. 1. The neutron scattering intensity in superfluid 4 He at 
saturated vapor pressure and T=1.3 K on a logarithmic scale 
as a function of wave vector Q and energy, from Ref . [l] . The 
line shows the dispersion of the quasiparticles, corresponding 
to phonons at small wave vectors Q and rotons near the local 
minimum at Qr ~ 1-93 A -1 and energy A/fcs ~ 8.6 K. 
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FIG. 2. Spin-echo scans as a function of coil displacement 
and temperature. With increasing temperature, the neutron 
polarization decreases (shown by the projected grey shadow) 
and the phase of the oscillation shifts. The statistical errors 
in the neutron count rates are smaller than the symbol size. 
The lines are fits by a cosine function. 



ror guide and a PG(002) monochromator, respectively, 
with the higher-order neutrons eliminated by a veloc- 
ity selector. Analysis of the final energy (6 meV) and 
polarization of the scattered neutrons were done with a 
PG(002) crystal assembly and a supermirror bender, re- 
spectively. The radio-frequency flippers were inserted in 
zero-field devices mounted before and after the sample. 
TRISP was operated in both boot-strap and non-boot- 
strap mode [5J in — + configuration and radio frequencies 
between 60 and 292.62 kHz, which cover spin-echo times 
15< r <190 ps. The superfluid 4 He sample, held in an 
aluminium cylinder of inner diameter 26 mm and volume 
30 cm 3 equipped with Cd spacers to reduce multiple scat- 
tering, was cooled to below 0.5 K using a 3 He insert in a 
pulse-tube cryocooler. Typical count rates for the roton 
was 60 counts per minute. 

The intensity of the scattered neutrons is recorded as 
a function of the length of one of the zero-field precession 
paths by scanning the displacement of one resonance coil. 
The oscillatory pattern can be fitted by a cosine function 
(see Fig. [2]) where the amplitude of the oscillatory com- 
ponent is a measure of the neutron polarization and gives 
the lifetime (inverse linewidth) of the excitation while the 
phase of the oscillations (position of maximum intensity) 
gives the relative energy of the excitation. 

The roton energy shift 



S(T) = A(T) - A(0) 



(1) 



was obtained in this way for temperatures between 0.5 
and 1.6 K using spin-echo times t of 50 and 115 ps. Since 
our measurements determine changes in the roton energy, 
and not the energy in absolute units, we take S(T)=Q at 
the lowest temperature of the measurements, T=0.5 K. 
To measure the roton linewidth T(T) (half width at half 
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FIG. 3. Temperature dependence of (a) the roton energy shift 
5(T) on a linear scale and (b) of the roton linewidth T(T) (full 
width at half maximum) on a logarithmic scale. Symbols are 
experimental data with one-si gm a erro r ba rs, solid line the 
best least-squares fit of Eqs. (12 1 and (13 1, dotted line the 



roton contribution, and dashed line the phonon contribution. 



maximum), we first verified that the decay of the ro- 
ton is exponential, by measuring the neutron polarization 
P(t,T) = Po exp[— r(T)r] for spin-echo times r between 
16 and 190 ps at six temperatures between 0.5 and 1.5 
K. Once the exponential decay was established, further 
measurements of P(r, T) at intermediate temperatures 
were performed for spin-echo times of 50, 115, and 190 
ps. The linewidth obtained from these data sets was 
put on an absolute scale using P obtained from direct 
beam calibration measurements and taking into account 
the depolarization due to beam divergence and dispersion 
curvature using the formalism of Ref. [H]. At T=0.5 K, 
where no measurable broadening of the roton line is ex- 
pected, we find a negligible linewidth (<2 mK), which 
confirms this calibration procedure. 

The experimental results are shown in Fig. [3] The 
roton linewidth, Fig. [3](b) , follows the temperature de- 
pendence given by roton-roton scattering that has been 
observed previously [3J HJ and agrees with the LK 
theory O [7] . On the other hand, the roton energy ini- 
tially increases with temperature, see Fig.^a), before it 
begins a monotonic decrease at T > 1 K. The change of 
sign of the roton energy shift provides unambiguous evi- 
dence that more than one interaction process contributes 
to it. In one earlier report lOj, the roton energy was 
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FIG. 4. Self-energy diagrams due to quasiparticle interac- 
tions in superfluid 4 He. Solid lines correspond to rotons and 
dashed lines to phonons. Red lines indicate thermally excited 
quasiparticles. Triangles and squares correspond to cubic and 
quartic vertices, respectively. Four-particle processes: (a) A 
roton is scattered by a thermal roton or (b) by a thermal 
phonon. In (c), the roton interacts with a uniform density 
of thermal phonons (the Hartree term). Three- particle pro- 
cesses: (d) decay and (e) coalescence of a roton with a thermal 
phonon. 



found to change more slowly than predicted by the LK 
theory [7 r However, no change of the sign of the energy 
shift was observed and no explanation of the behavior 
was found at that time. 

In order to explain these experimental findings, we use 
a diagrammatic approach [12] to evaluate the contribu- 
tions from the lowest-order self-energy diagrams shown in 
Fig. [4j In contrast to previous work, we explicitly take 
into account the interaction of a roton with thermally 
excited low-energy phonons, and include both three- and 
four-particle processes. 

The 4PP interaction in the diagrams of Figs.[5]Ja)-[4][c) 
is given by 



V 4 



(2) 



with p + q = p' + q'. The quartic vertex Vp q;p < q ' is a 
constant for roton-roton scattering and U pq;P ' q ' ~ \fq~q 1 
for roton-phonon interaction processes jH [7j, where q 
and q' are the phonon momenta. We assume that the 
energy of quasiparticles with momentum p close to the 
roton minimum pr can be parameterized by the Lan- 
dau spectrum e p = A + (p — pr) 2 / (2m*) with m* being 
the roton mass, while quasiparticles in the phonon region 
with momentum q are assumed to have a linear dispersion 
uj q = cq. We use numerical values of A/fee = 8.622 K 
PQ, Pr/K = Qr = 1-931 A" 1 p, n = m*/mn e = 0.144 
PQ, c = 238.3 m/s p3], and p = 0.14513 g/cm 3 [M], and 
express linewidths and energy shifts in units of Kelvin. 

Roton-roton scattering, i.e., the scattering of a roton 
by a thermally excited roton, is depicted as a self-energy 
diagram in Fig. Qa). The roton energy shift and the 
linewidth due to this process using Eq. ([2| are both pro- 



portional to the number of thermally excited rotons, 
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as first shown by LK [2]. The minor correction factor in 
brackets is due to small deviations of the roton dispersion 
from a parabolic form, with a^JpZ — 0.0603 K -1 / 2 j71ll5j. 

Roton-phonon scattering, i.e., the scattering of a roton 
by a thermally excited phonon, is represented by the di- 
agram in Fig. |4jb), which is essentially the same as that 
for roton-roton scattering. However, since phonons are 
gapless, the linewidth and the energy shift of a roton due 
to this scattering process depend on temperature via a 
power law. The linewidth using Eq. ^ is 

r p (T) = |Vp q;pV | 2 n q (l + ?v +7 V) 
q,q' 



-J <5(e p +w q -w q / -e p /), (4) 



which yields the low-T asymptote 
T P (T) 
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where A is a dimensionless constant characterizing the 
roton-phonon interaction. It can be evaluated in terms 
of density fluctuations of the phonon field in the quantum 
hydrodynamic limit as [5] 
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The numerical value of A is not known with precision. 
Early estimates indicated a negative value —0.1 [16], 
while more recent neutron scattering measurements give 
values in the range -0.38 to +1.35 [10 [IS]. The un- 
certainty comes mostly from the roton mass and from 
the density dependence of the roton energy A(p). Using 
values of A in this range, the prefactor of the T 7 -term 
in Eq. ^ turns out to be so small that this process 
is negligible compared to that of the roton-roton scat- 
tering [Eq. (|3|] down to T — 0.5 K, where both mech- 
anisms yield an undetectable linewidth of the order of 
10~ 6 K. Qualitatively, the interaction of the massless 
long-wavelength phonon with a massive almost immo- 
bile roton is analogous to Rayleigh scattering. Multi- 
plying the corresponding cross-section a oc a; 4 ~ T 4 by 
the density of thermally excited phonons, N p h{T) cx T 3 , 
leads to a T 7 -law. The roton energy shift evaluated from 
the roton-phonon diagram in Fig. |4jb) has an equivalent 
T 7 contribution, which can also be neglected compared 
to the T 4 contributions from other interaction processes 
that are discussed next. We also note that a T 7 power 
law for the energy shift is not consistent with our exper- 
imental data. 

The Hartree term, a 4PP illustrated in the diagram of 
Fig. |4^c), contributes only to the real part of the self- 
energy. It corresponds to an energy renormalization of a 
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roton propagating in a gas of thermally excited phonons. 
Straightforward calculation of this diagram using Eq. §} 
gives the energy shift 



a 4 (r) = £ ff (p) « £V P . 



with a low-temperature asymptote of the form 
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Three-particle processes are depicted in Figs. (4^d) and 
|4|e), where the first diagram corresponds to decay of a 
roton into a roton and a phonon and the second diagram 
describes absorption of a thermal phonon by a roton. 
For the dispersion shown in Fig. [T] these three-particle 
processes do not simultaneously conserve energy and mo- 
mentum and do not contribute to the roton lifetime, i.e., 
the resulting roton self-energies have no imaginary part. 
However, the real part of the self-energy, which corre- 
sponds to virtual transitions into higher-energy states, is 
non-zero and gives rise to a temperature-dependent en- 
ergy shift. The 3PP interaction is given by 



V, = 



E 

p^q.p' 



t t 



(9) 



with p = q + p'. For interaction processes involving one 
phonon q, the cubic vertex scales as Vp ;q p< cx ^fq (for the 
full expression, see Ref. [5]). The net contribution of the 
3PP diagrams in Figs. Qd) and[4|e) to the temperature 
dependent part of the roton energy is 
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At low-temperatures, Eq. (10) gives a T law 
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(10) 



(11) 



similar to the Hartree term of Eq. ([8]) but with a differ- 
ent prefactor. While the contribution of Eq. (11) to the 
energy shift is positive, that of Eq. Q depends on the 
sign of the roton-phonon interaction. 

We conclude from these theoretical considerations that 
the dominating process contributing to the temperature- 
dependent linewidth of the roton, T(T), is the roton- 
roton scattering, while the shift of the roton energy, 5(T), 
is due to both roton-roton scattering and the renormal- 
ization of the roton energy due to the background of ther- 
mally excited phonons. Thus, we obtain 

T(T) = 7P T 7 + lR VT \l+ay/0\ e~ A ( T >/ T (12) 



and 



S(T) = S P T 4 - S R Vf \l + a^/J/r] e 



-A(T)/T 



(13) 



The four parameters 7p, jr, 6p, and Sr were de- 
termined by fitting Eqs. ( fl2| and (13) self-consistently 
[since A(T) depends on Eq. (13)] to our experimental 
data. The lines in Fig. [3] show the resulting least-squares 
fit. Clearly, the T 4 -term explains the upturn of the ro- 
ton energy at the lowest temperatures, before the down- 
turn above 1 K, where roton-roton scattering dominates. 
For the roton linewidth, inclusion of the T 1 term due 
to roton-phonon scattering processes does not improve 
the description of the experimental data in the measured 
temperature regime, T > 0.5 K. 

The best fit gives roton-roton scattering contributions 
of 7 _r = 49.6 ± 1.1 K 1 / 2 and S R = 25.6 ± 3.0 K 1 / 2 , 
in a good agreement with previous experimental deter- 
minations, which are in the ranges jr « 42 — 47 and 
S R w 19 - 25, respectively \M EE] • The T 7 roton- 
phonon contribution to the roton linewidth is not signif- 
icant, as shown in Fig. [3^b) , where we have plotted as a 
dashed line an estimate of this contribution, taking the 
largest value of the roton-phonon interaction parameter 
compatible with the measured density dependence of the 
roton parameters, A = 1.35. On the other hand, the 
T 4 -contribution to the roton energy is clearly significant 
in the experimental fit, see Fig. [3]Ja), yielding a value of 
8p = 8.0± 1.5 x 10 -3 K -3 . Our theoretical estimate from 
Eqs. (|| and (n} with A = 1.35 gives 8$ = 5x 10~ 3 K" 3 , 
close to the experimentally observed value. Vertex cor- 
rections to the diagrams in Figs.[4|d) and[4je), which are 
beyond the scope of the present theoretical calculations, 
increases the value of 63 and would further improve the 
agreement. 

In conclusion, our work shows that at sufficiently low 
temperatures, the roton energy is affected not only by 
scattering from thermally excited rotons but also by the 
presence of a "thermal bath" of excited phonons. The ef- 
fect is particularly visible since the two processes give rise 
to different signs of the energy shift. The magnitude of 
the initial upward shift of the roton energy suggests that 
it has contributions from both the three-particle cubic 
vertex terms of Figs.Qd) and[4je) and from the quartic 
Hartree-term of Fig. Hf c) . We further conclude that the 
roton-phonon interaction parameter A is likely to be pos- 
itive, and not negative as was assumed until now. This 
implies that the effective roton-phonon interaction is re- 
pulsive, in contrast to the roton-roton interaction, which 
is attractive |19j . The roton-phonon interaction does not 
give any observable roton linewidth. The present results 
show the importance of both cubic and quartic vertices 
in superfluid 4 He, and we expect that this will have a 
wide-ranging impact for Bose-condensed systems in gen- 
eral, covering areas such as excitons in semiconductors, 
magnon decay [20H22) . non-linear effects in spintronics 
[23l [24] . and strongly interacting ultracold Bose gases 
[25]. 
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